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production at hadron colliders, we present here the calculations of gluon jet 
Q\ ■ production at hadron colliders in the two-gluon exchange parameterization of 

the Pomeron model. We use the helicity amplitude method to calculate the 
cross section formula. We find that for the gluon jet production the diffractive 
Q_i! process is related to the differential off-diagonal gluon distribution function in 

the proton. We estimate the production rate for this process at the Fermilab 
Tevatron by approximating the off-diagonal gluon distribution function by 
the usual diagonal gluon distribution. 
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I. INTRODUCTION 



In recent years, there has been a renaissance of interest in diffractive scattering. These 
diffractive processes are described by the Regge theory in terms of the Pomeron (IP) exchange 
The Pomeron carries quantum numbers of the vacuum, so it is a colorless entity in QCD 
language, which may lead to the "rapidity gap" events in experiments. However, the nature 
of Pomeron and its reaction with hadrons remain a mystery. For a long time it had been 
understood that the dynamics of the "soft pomeron" is deeply tied to confinement. However, 
it has been realized now that how much can be learned about QCD from the wide variety 
of small- x and hard diffractive processes, which are now under study experimentally. In 
Refs. the diffractive J/ip and T production cross section have been formulated in 

photoproduction processes and in DIS processes in perturbative QCD. In the framework 
of perturbative QCD the Pomeron is represented by a pair of gluon in the color-singlet 
sate. This two-gluon exchange model can successfully describe the experimental results 
from HERA §. 

On the other hand, as we know that there exist nonfactorization effects in the hard 
diffractive processes at hadron colliders First, there is the so-called spectator effect 

||, which can change the probability of the diffractive hadron emerging from collisions intact. 
Practically, a suppression factor (or survive factor) u Sp n is used to describe this effect [ID 



Obviously, this suppression factor can not be calculated in perturbative QCD, which is now 
viewed as a nonperturbative parameter. Typically, the suppression factor Sp is determined 
to be about 0.1 at the energy scale of the Fermilab Tevatron |J. Another nonfactorization 
effect discussed in literature is associated with the coherent diffractive processes at hadron 
colliders 0, in which the whole Pomeron is induced in the hard scattering. It is proved in 
that the existence of the leading twist coherent diffractive processes is associated with a 
breakdown of the QCD factorization theorem. 

Based on the success of the two-gluon exchange parametrization of the Pomeron model 
in the description of the diffractive photoproduction processes at ep colliders J2],|3|,|5j , we may 
extend the applications of this model to calculate the diffractive processes at hadron col- 
liders in perturbative QCD. Under this context, the Pomeron represented by a color-singlet 
two-gluon system emits from one hadron and interacts with another hadron in hard process, 
in which the two gluons are both involved (as shown in Fig. 1). Therefore, these processes 
calculated in the two-gluon exchange model are just belong to the coherent diffractive pro- 
cesses in hadron collisions. Another important feature of the calculations of the diffractive 
processes in this model recently demonstrated is the sensitivity to the off-diagonal parton 



distribution function in the proton 



Using this two-gluon exchange model, we have calculated the diffractive J/ip production 
12fl , quark jet production massive muon pair and W boson productions [fl5[] , and 



direct photon production |]TB[ in hadron collisions. In this paper, we will further calculate 
the gluon jet production at large transverse momentum in the coherent diffractive processes 
at hadron colliders by using the two-gluon exchange model. In the calculations of Refs. 



I2] , p!3|jr5| , there always is a large mass scale associated with the production process. That 
is Mip for J/ip production, m c for the charm jet production, M 2 for the massive muon 
production (M 2 is the invariant mass of the muon pair) and for W boson production. 
However, in the gluon jet production process as well as the light quark jet production process, 
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there is no large mass scale. So, for these processes, the large transverse momentum is needed 



to guarantee the application of the perturbative QCD. Furthermore, in [|14| we show that 
the light quark jet production in the two-gluon exchange model has a distinctive feature 
that there is no contribution from the small l\ region {l\ < k^) in the integration of the 
amplitude over l\. (The similar behavior has also been found for the diffractive light quark 
photoproduction process |17|.) So, the expansion (in terms of l\jM\) method used in Refs. 
P^ , P^B|JT5| can not be applied to the calculations of gluon jet production. In the following 
calculations, we will employ the helicity amplitude method to calculate the amplitude of the 
diffractive gluon jet production in hadron collisions. We will show that the production cross 
section is related to the differential (off-diagonal) gluon distribution function in the proton 
as that in the diffractive light quark jet production process ||14|| . (On the other hand, we 
note that the cross sections of the processes calculated in Refs. [|1^,[13|,|T^] are related to the 
integrated gluon distribution function in the proton). 

Diffractive gluon jet production can come from two types partonic processes: one is the 
quark initiated process (Fig. 2), and the other is the gluon initiated process (Fig.3). 

The diffractive production of heavy quark jet at hadron colliders has also been studied by 
using the two-gluon exchange model in Ref. ||18|| . However, their calculation method is very 
different from ours []. In their calculations, they separated their diagrams into two parts, 
and called one part the coherent diffractive contribution to the heavy quark production. 
However, this separation can not guarantee the gauge invariance 0]. In our approach, 
we follow the definition of Ref. 0, i.e., we call the process in which the whole Pomeron 
participants in the hard scattering process as the coherent diffractive process. Under this 
definition, all of the diagrams plotted in Fig. 2 and Fig.3 for the partonic processes contribute 
to the coherent diffractive production. 

The rest of the paper is organized as follows. In Sec. II, we will give the cross section 
formula for the partonic process in the leading order of perturbative QCD. In this section 
we employ the helicity amplitude method to calculate two partonic processes, qp — > qgp 
and gp — > ggp. In Sec. Ill, we estimate the production rate of diffractive gluon jet at the 
Fermilab Tevatron by approximating the off-diagonal gluon distribution function by the 
usual diagonal gluon distribution function in the proton. We also compare the contributions 
from different partonic processes to the diffractive dijet production at the Tevatron. And 
the conclusions will be given in Sec. IV. 

II. THE CROSS SECTION FORMULA FOR THE PARTONIC PROCESS 

A. qp — ► qgp process 

For the partonic process qp —>■ qgp, in the leading order of perturbative QCD, there are 
nine diagrams shown in Fig. 2. The two-gluon system coupled to the proton (antiproton) in 
Fig. 2 is in a color-singlet state, which characterizes the diffractive processes in perturbative 
QCD. Due to the positive signature of these diagrams (color- singlet exchange), we know 



x For detailed discussions and comments, please see [13] 



3 



that the real part of the amplitude cancels out in the leading logarithmic approximation. 
To get the imaginary part of the amplitude, we must calculate the discontinuity represented 
by the crosses in each diagram of Fig. 2. 

The first four diagrams of Fig. 2 are the same as those calculated in the diffractive direct 
photon production process at hadron colliders But, due to the existence of gluon- 

gluon interaction vertex in QCD, in the partonic process qp — > qgp, there are additional five 
diagrams (Fig.2(5)-(9)). These five diagrams are needed for a complete calculation in this 
order of QCD. 

In our calculations, we express the formulas in terms of the Sudakov variables. That is, 
every four-momenta fcj are decomposed as, 

ki = a { q + (3ip + k iT , (1) 

where q and p are the momenta of the incident quark and the proton, q 2 = 0, p 2 = 0, and 
2p ■ q = W 2 = s. Here s is the cm. energy of the quark-proton system, i.e., the invariant 
mass of the partonic process qp — > qgp. on and $ are the momentum fractions of q and p 
respectively. kiT is the transverse momentum, which satisfies 

k iT ■ q = 0, k iT ■ P = 0. (2) 

All of the Sudakov variables for every momentum are determined by using the on-shell 
conditions of the momenta represented by the external lines and the crossed lines in the 
diagram. The calculations of these Sudakov variables are similar to those in the diffractive 
light quark jet production process gp — > qqp fl4|, and we can get the Sudakov variables 



of every momentum for the process qp — > qgp from the relevant formulas of |L4|. In the 



following, we list all of the Sudakov variables for the diffractive process qp — > qgp. 
For the momentum u, we have 

Ml 2 , , 

oi u = 0, p u = xjp = — — , u T = t = 0, (3) 

where M\ is the invariant mass squared of the diffractive final state including the light 
quark and antiquark jets. For the high energy diffractive process, we know that M 2 <C s, so 
we have (3 U (xjp) as a small parameter. For the momentum k, 

k 2 

a k (l + a k ) = -tJ^-> Pk = -akPu, (4) 

where k? is the transverse momentum of the out going quark jet. For the loop mentum I, 
because the results for Pi are not the same for the nine diagrams of Fig.2, we get its value 
from the formula of Ref. [14| for the relevant diagram. The results are 

Oil = , 

s 

2(kT,h)-^ for m ! 2 6 
a k s 



2(k T ,l T )+l 2 T 
(l + a k )s 
M 2 x -l 2 T 



for Diag.5, 7, 8, 

for Diag.3, 4, 9, (5) 
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where (h?, It) is the 2- dimensional product of the transverse vectors hr and It- 

Using these Sudakov variables, we can give the cross section formula for the partonic 
process qp — > qgp as, 

da(qp^qgp) dM\d 2 k T da k fcf, 2 

It = IG^lG^Mi 6{M1 + Qk) + Ml )i: 11 ' (6) 

where A is the amplitude of the process qp — > qgp. We know that the real part of the 
amplitude A is zero, and the imaginary part of the amplitude A(qp — > qgp) for each diagram 
of Fig. 2 has the following general form, 

ImA = C F (T«) J x Ui{u - fc)lWg), (J) 

where Cp is the color factor for each diagram, a is the color index of the incident gluon. T M 
represents some 7 matrices including one propagator. F in the integral is defined as 

F = ^-g 3 J(x',x";l T ), (8) 

where 

f/V t"-1 2 \ - d G ( x '> x "i 1 t) / q x 



where the function G(x', x"\ k T ) is the so-called off-diagonal gluon distribution function HIT 
Here, x' and x" are the momentum fractions of the proton carried by the two gluons. It is 
expected that at small x, there is no big difference between the off-diagonal and the usual 



diagonal gluon densities ||19|| . So, in the following calculations, we estimate the production 
rate by approximating the off-diagonal gluon density by the usual diagonal gluon density, 
G(x',x";Q 2 ) « xg(x,Q 2 ), where x = x F = M\js. 

In Hl6| , we calculate the cross section for the diffractive photon production process qp — > 
'yqp by directly squaring the partonic process amplitude. However, in the calculations here 
for the partonic process qp —* qgp because there are additional five diagrams contribution, 



it is not convenient to directly square the amplitude. Following Ref. [[[J], we calculate the 
amplitude by employing the helicity amplitude method pO|j2~If . Furthermore, we will show 
that by using the helicity amplitude method we can reproduce the cross section formula for 



the diffractive photon production process [16 



For the massless quark spinors, we define 

u ± (p) = ^=(l± l5 )u(p). (10) 

For the polarization vector of the outgoing gluon (its momentum is k + q), following the 
method of Ref. 1 12111, we find that it is convenient to choose 



= N e \^+ i) i T 75)+ i> m+ m ± 75)]. (ii) 

The normalization factor N e equals to 
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N e = — _. (12) 



s^2k 2 T 

With this definition ([0]), we can easily get the scalar products between the four-momenta 
and the polarization vector e as 

k 2 s 

e • p — 0, e ■ q = N e — - — , e-kT = —N e k T s, e ■ It = —N e (hr,lT)s. (13) 
1 + oi k 

The helicity amplitudes for the processes in which the polarized Dirac particles are 
involved have the following general forms 



u±{p f )Qu±{pi) = (14) 

where n is an arbitrary massless 4-vector, which is set to be n = p in the following calcu- 
lations. Using this formula ([14]) , the calculations of the helicity amplitude *A(Ai, A2, A3) for 
the diffractive process qp — > qgp is straightforward. Here Ai represents the helicity of the 
incident quark; A 2 and A 3 represent the helicities of the outgoing gluon and quark respec- 
tively. In our calculations, we only take the leading order contributions, and neglect the 

M 2 

higher order contributions which are proportional to p u = — p because in the high energy 
diffractive processes we have ^<1. 

For the first four diagrams, to sum up together, the imaginary part of the amplitude 
A{+,+,+) is 

d% f . , „ 2 2 -lk T ~{l + a k ){k T ,l T ). 



lmA 12U (+,+,+) = al(l + a k )Afx J —Lf{ x \ x "-l 2 T ){ 



(/2)2- ' ' ^9 36 {k T - (1 + ak )l T )2 - 

(15) 



where | and are the color factors for Diags.1,4 and Diags.2,3 respectively, and Af is 
defined as 

A/~= , 35 „^V;- (16) 



'— 1a k k T 

The other helicity amplitudes for the first four diagrams have the similar forms as (|i~5|) , 
W 1234 (-, -, -) = W 1234 (+, +, +), 

W 1234 ( + , -, +) = W 1234 (-, +, -) = Zll na 1234 (+) + +) _ (17) 

ttfc 

These amplitude expressions Eq. ( |T5| ) can also serve as the calculations of the amplitude 
for the diffractive direct photon production process qp — > q'-fp [[3IJ except the difference on 
the color factors.0 In the direct photon process, the color factors for these four diagrams 



2 In Ref. p6|, we did not employ the helicity amplitude method. If we use the amplitude expressions 



Eqs.(15) and ( JT7| ) (correct the color factors) to calculate the photon production process qp — ► qjp, 
we can get the same result as that in filf . This can be viewed as a cross check for the methods we 
used in the calculations. 
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are the same (they are all |). It is instructive to see what is the consequence of this 
difference. We know that the amplitude of the diffractive process in Eq. (0) must be zero 
in the limit l\ — > 0. Otherwise, this will lead to a linear singularity when we perform the 
integration of the amplitude over l\ due to existence of the factor \j{\\) 2 in the integral of 
Eq. (0) ||13|| . This linear singularity is not proper in QCD calculations. So, we must first 



exam the amplitude behavior under the limit of l\ — > for all the diffractive processes in 
the calculations using the two-gluon exchange model. From Eq. (|T^), we can see that the 
amplitude for the diffractive direct photon production process qp — > q^p is exact zero at 
l\ — > 0. However, for the process qp — > qgp the amplitude for the first four diagrams is not 
exact zero in the limit l\ — > due to the inequality of the color factors between them. So, 
for this process there must be other diagrams in this order of perturbative QCD calculation 
to cancel out the linear singularity which rises from the first four diagrams. The last five 
diagrams of Fig.2 are just for this purpose. 

Finally, by adding up all of the nine diagrams of Fig.2, the imaginary parts of the 
amplitudes are 



a 2 



Ira4(+, +, +) = IrrU(- -, -) = -jM x T, 
InU(+, -, +) = InU(- +, -) = -^N x T, (18) 



where 



t f rf2/T t( ' " ^ J 1 +ak) 2 {k T J T ) - (1 + ak)l 2 T , (fcr, h) + l\ 
1 = / /,o no T\x , x ; l T )\ (1 + afc) — =; 

(k,T, It)— It , i{^t-, It) — otkl T -, 

—dk — h Oi k — =; =; J- (19) 

(hr ~ It) 2 ' (k T - akh) 2 

From the above results, we can see that in the integration of the amplitude the linear 
singularity from different diagrams are canceled out by each other, which will guarantee 
there is no linear singularity in the total sum. 

Another feature of the above results for the amplitudes is the relation to the differential 
off-diagonal gluon distribution function f{x',x";l^). However, as mentioned above that 
there is no big difference between the off-diagonal gluon distribution function and the usual 
gluon distribution at small x, so we can simplify the integration of ( |19|) by approximating 
the differential off-diagonal gluon distribution function f(x', x"\ I 2 -) by the usual diagonal 
differential gluon distribution function f g (x; l?)- 

After integrating over the azimuth angle of It, the integration T will then be 

dlj> 2 1 + a k 1 kj, — (1 + ak)l T ,1 k\ — l\ 



t = * i£rM* - n T 2I ) + (5 - 



(llY J9K ' in 9 v 2 2\k 2 T - {I + a k )l 2 T \' v 2 2\k 2 T -l 2 T 

1 k\ — a k l 2 

2 2 | Otfclrp 



2 , £, )]• (2°) 



In the above integration, if I 2 < k T /(l + a k ) 2 the first term of the integration over l\ will 
be zero; if If < k\ the second term will be zero; if l\ < k^/aj. the third term will be zero. 
So, the dominant regions contributing to the three integration terms are I 2 ~ k\j{\ + a k ) 2 , 
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If ~ fcy, and If ~ k^/al respectively. Approximately, by ignoring some evolution effects 
of the differential gluon distribution function f g {x]lj) in the above dominant integration 
regions, we get the following results for the integration T, 

T = £[/,(x; 4) + ^/.ta + oj/.fe |)]. (21) 

Obtained the formula for the integration T, the amplitude squared for the partonic 
process qp — > qgp will be reduced to, after averaging over the spin and color degrees of 
freedom, 

lAl ~ 24 Ml{l + a k ) S 111 ■ {22) 
And the cross section for the partonic process qp — > qgp is 



da(qp -> qgp) 

\t=o 



dt JM^>4k T 



/ dMxdk"^dak[5(ak — «i) + 5(ak — a 2 )] 



«2 1 in 2 , (23) 



96(M|-) 2 1 + a fc 



where ai )2 are the solutions of the following equations, 



k. 



a(l + a) + = 0. (24) 

The integral bound M x > Ak\ in ( |23"D shows that the dominant contribution of the inte- 
gration over M\ comes from the region of M\ ~ Ak 2 -.. Using Eq. (^), this indicates that 
in this dominant region is of order of 1. So, in the integration T the differential gluon 
distribution function f g (x; Q 2 ) of the three terms can approximately take their values at the 
same scale of Q 2 = k\. That is, the integration T is then simplified to 

T = 9p-/»fa ^X 1 + a *)( 10 - 7c ^ + 10 «fc)- ( 25 ) 

Numerical calculations show that there is little difference between the cross sections by 
using these two different parametrizations of T, Eq. ( pl|) and Eq. (^). So, in Sec. IV, we 
use Eqs. (1231) and (E5f) to estimate the diffractive production rate at the Fermilab Tevatron. 



B. gp — > ggp process 

For the partonic process gp — > ggp, there are twelve diagrams in the leading order 
contributions as shown in Fig.3. The first nine diagrams are due to the existence of the 
three-gluon interaction vertex, and the last three diagrams are due to the existence of the 
four-gluon interaction vertex. But it will be shown in the following calculations, the last 
three diagrams do not contribute under some choice of the polarizations of the three external 
gluons. 



S 



The Sudakov variables can be calculated by the similar method used in the last sub- 
section. And those Sudakov variables of the momenta u, k for the process gp — > ggp are 
the same as those in the last subsection. For the loop momentum I, the relevant Sudakov 
variables for each diagram are 

®l = , 

S 

Pi = gfeM^gl , for Diag.l, 4, 6, 10, 



2(k T ,l T ) + l 2 T 
(1 + a k )s 
M\-l\ 



for Diag.2, 3, 5, 11, 

for Diag.7, 8, 9, 12. (26) 



And also, we can express the cross section formula for the partonic process gp — > ggp in 
the following form, 

dv(gp^ggp), dM\d 2 k T da k 2 

It lt=0 = m«/ W1 + ak) + 11 ' (27) 

where A is the amplitude of the process gp — > ggp. We know that the real part of the 
amplitude A is zero, and the imaginary part of the amplitude A(gp — > ggp) for each diagram 
of Fig. 3 has the following general form, 

ImA = C F f abc I ^G(ei, e 2 , e 3 ) x F, (28) 
J {l T ) 

where Cp is the color factor for each diagram, a, b, c are the color indexes for the incident 
gluon and the two outgoing gluons respectively, and f a b c are the antisymmetric SU(3) struc- 
ture constants. G(ei, e 2 , 63) represents the interaction part including one propagator for the 
first nine diagrams, where e\, e 2 , are the polarization vectors for the incident gluon and 
the two outgoing gluons. F in the integral is the same as that in Eq.(§). 
The color factors Cp for the twelve diagrams are 



C F = 


1 

2' 


for Diag 




7, 


C F = 


1 

~2' 


for Diaj 


5-2, 




C F = 


1 

~4' 


for Diaj 


r.3, 


6, 9, 


C F = 


1 

4' 


for Diag 


;-4, 


5, 8, 


C F = 


3 
4 : 


for Diag 


;.10 


, 11, 12 



(29) 

Following the calculation method used in the last subsection, we employ the helicity 
amplitude method to calculate the amplitude Eq.(^). For the polarization vector of the 
incident gluon, which is transversely polarized, we choose, 
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^ = -^(0,1,^,0). (30) 



For the two outgoing gluons, we choose their polarization vectors as |2T[ 

A ±] =N e U+ f) ^(i T75 )+^(# + 0(i± 75 )], 

A ±] = N e [{i- ft 4 i>{\ T 75 )+ i> 4(i- W(i ± 75)]. (3i) 



The normalization factor N e has the same form as in Eq.(O). Under the above choice of 
the polarization vectors for the external gluons, we can easily find that they are satisfied the 
following equations, 

P ■ ei = p ■ e 2 = p ■ e 3 = 0. (32) 

With these relations, we can further find that the last three diagrams do not contribute to 
the partonic process gp — > ggp. 

For the first nine diagrams, there are two helicity amplitudes among the eight helicity 
amplitudes do not contribute in the context of the above choice of the polarizations of the 
external gluons, i.e., 

lmA{+, +, +) = IrrU(- -, -) = 0. (33) 

In the expression of the amplitude A(X(ei), Xfa), A(es)), A denote the helicities for the three 
gluons respectively. The other six helicity amplitudes are divided into the following three 
different sets, 

lmA{+, -, -) ~ lmA{-, +, +), 
Inu4(+, -, +) ~ lmA{-, +, -), 

lmA{+,+,-) ~InU(--+). (34) 

For the first helicity amplitudes set, Ihl4.(±, qp, T)> to sum up all of the nine diagrams, 
we get 

InU(±, =F, T) = A/"Vef } • k T l, (35) 



where M' is defined as 



And the integration X is 



3 s 



= 7T2 9sfa b c (36) 



4fc 



2 i> s 
T 



1 r d 2 l T 2 k T + (k T ,l T ) k\ - (k T , It) 

T<W S(X • X ' ' T>H1 + ^ (kr + hY + ^ (kr-T T r- 



71 



+ (1 + ak) 2k T -(l + a k )(k Tj l T ) _ a2 k T -a k (k Tj l T ) 

(k T - (1 + a k )l T f {k T -a k l T f 
1 f d 2 l T ,,. 1 2^u^ , Ak T ,l T ) + l T (k T ,l T )-l T 



» J W? fiX ' X ' ' T)[(1 + ^ Ck T + l T y + (kr - It, 



2 



2 (k T ,l T ) - (l + a k )l T 2 (k T ,l T ) - a k l\ 
a k) — — 7TT^ + a k— FT^J- ( 38 ) 



(kr - (1 + OL k )l T ) 2 " (k T - oc k l T ] 
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From the above equations, we can check that there is no linear singularity at the limit of 
l T — > in the integration of the amplitude over the loop momentum. The first term of 
the integration X in Eq.([J7|) comes from the contribution of Diag.3; the second term comes 
from Diag.4; the third term comes from Diag.6 and Diag.9; the last term comes from Diag.5 
and Diag.8. The contributions from Diags.l, 2 and 7 are canceled out by each other. From 
fl37p, we can see that the linear singularities coming from the four terms are canceled out 
by each other. The final result for the amplitude is now free of linear singularity. We must 
emphasize here that only the total sum of the contributions from all of the diagrams is free 
of linear singularity. The separation of these diagrams will cause linear singularity. 

Following the argument in the last subsection of the calculation for the partonic pro- 
cess qp — > qgp, we can approximate the differential off-diagonal gluon distribution function 
f(x', x"; l T ) by the usual diagonal differential gluon distribution function f g (x; l T ) to further 
simplify the integration of X. After integrating over the azimuth angle of It, this integration 
will then be 



1 = / wki'to *m - - (1 + aM 



(llY J9y i inK 2 2\k T -l T \' K K>y 2 2\k 2 T - (l + a k )l T \ 



1 k\ — a k l 2 

2 2 | krp Qt fa J,rp 



+»*(n- o.S ( 39 ) 



The above equation shows that the integration X here has the similar behavior as that of the 
integration T of Eq.(|20|) in the last subsection. So, the three terms of the above integration 
X are dominantly contributed from the integral regions of l T as l\ ~ k T /(l + a k ) 2 , l\ ~ k T , 
and I 2 ~ k T /a\ respectively. Approximately, we may also ignore the evolution effects of the 
differential gluon distribution function f g (x;l T ) in the above dominant integration regions, 
and so the integration X is reduced to 

I = 4) " (1 + u k ff g (x- I ' 4 2 ) + a\f g {x- §)]. (40) 

lirp yl_ -\- Cx k j CXy. 

For the second helicity amplitudes set, Inu4.(±, =p, ±), the calculations are more compli- 
cated, and the contribution from Diag.3 is 



W 3 (±, T , ±) = + a k f J m^f(x', x"; l 2 T ) 



a k k T ^ ] ■ (k T + l T ) + (k\ + (hr, l T ))^ ] ■ k T 



(k T + i T ) 2 



(41) 



The contribution from Diag.4 is 



lmA\±,T,±) =M'a k {l + a k ) f ^ f ( x ' , x " ; l 2 T ) 

J [i T ) 

a k k T ^ ] ■ (fc r - l T ) + (k 2 T - (k T , / T ))ef } ■ k T 
{k T - h) 2 

The contributions from Diag.5 and Diag.8, to sum up together, are 



(42) 
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Inu4 58 (±, T ,±) = -M'a k (l + a k ) I ^f( x ',x";l 2 T ) 

J {l T ) 

a k k T e^ ] ■ (k T - a k l T ) + {k\ - a k (k T} h))^ ■ k T 
(k T - a k l T ) 2 

The contributions from Diag.6 and Diag.9, to sum up together, are 

Inu4 69 (±, T ,±) = N'(l + a k ) 2 I f^f(x',x";l T ) 

a k k 2 T ^ ] ■ (kr ~ (1 + ajS) + (4 - (1 + a k )(k T , h))^ ■ k T 

{k T -{l + a k )l T y 



(43) 



(44) 



The contributions from other three diagrams (Diag.l, Diag.2 and Diag.7) are canceled out 
by each other. From the above results Eqs.( |4l^]4*4] ), we can see that every term has linear 
singularity at the limit of l T — > in the integration of the amplitude over l T , while their 
total sum is free of the linear singularity. 

Following the procedure as we do for the helicity amplitude A(±, =F, =f) in the above, 
we can approximate the off-diagonal gluon distribution function f(x',x";l T ) by the usual 
diagonal differential gluon distribution function f g (x; l T )- After integrating over the azimuth 
angle of It, to sum up all of Eqs. (f41] - pi4|) , we get the helicity amplitude, 

W(±, =p, ±) = M'ir{\ + a k ) 2 e^ ] • k T X, (45) 



where I is the same as Eq.(^) and then Eq.(^D|) under the same approximation. 

For the third helicity amplitudes set, Iul4(±, ±, =)=), the calculations are similar to the 
calculations of Iul4(±, =)=, ±). The contribution from Diag.3 is 



Inu4 3 (±, ±, T ) = -M'a k {\ + a k ) J jj^fW, x"; l 2 T 



(Pin 

(l T J 

(1 + a k )k T et } • (kr + W) - {k 2 T + (k T , h))^ ■ k T 



2 



(46) 



The contribution from Diag.4 is 

d 2 h 



W 4 (±,±, T ) =M'a\ J ^- 2 f(x',x";l 2 T ) 



(1 TJ 

(1 + Ofcjgeg ■ (kr - l T ) - (k 2 T - (k T , l T ))^ ■ k T 
(k T - W) 2 

The contributions from Diag.5 and Diag.8, to sum up together, are 

Inu4 58 (±,±, T ) = -N'a 2 k f ^f(x',x";l T ) 

(1 + akjk^e^ ■ (k T - a k l T ) - (k\ - a k (k T , h))^i ■ k T 



(47) 



(k T - oi k l T ) 



2 



(48) 
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The contributions from Diag.6 and Diag.9, to sum up together, are 

lmA m (±,±,T) = M'a k {l + a k ) f f ( x ' , x » ; 

J {i T ) 

(1 + a k )k T ^ ■ (k T - (1 + a k )h) - (k T - (1 + a k )(k T , l T ))tf* ■ k 7 



(k T - (1 + a k )l 



T 



\2 



(49) 



And also, we find that the contributions from other three diagrams (Diag.l, Diag.2 and 
Diag.7) are canceled out by each other, and the total sum of Eqs. (f46|-fi9|) is free of the linear 
singularity. If we approximate the off-diagonal gluon distribution function f(x',x";l T ) by 
the usual diagonal differential gluon distribution function f g (x;l T ), and integrate over the 
azimuth angle of It, their sum will lead to a similar result as in Eq.( 



lmA{±, ±, =f) = Af'Tia 2 ^ ■ k T l. (50) 



By summing up all of the helicity amplitudes Eqs.(^), ([45]) and |5(J), we will get the 
amplitude squared for the partonic process gp — > ggp, after averaging over the spin and 
color degrees of freedom, 



-2 277T 2 «3(47T) 3 S 2 ^Jfi 



w = _ (51) 



And the cross section for the partonic process gp —> ggp is 



dajgp^ggp) = r 2 2 27oty _ % 23 1 

dt 1 A4>44 x T Z2M\k T M\ } 1 1 L _ 4k% 



Following the same argument in the last subsection for the calculations of the partonic 
process qp — > qgp, we see that the dominant contribution of the integration over M\ comes 
from the region of M\ ~ 4k T , where the differential gluon distribution function f g (x; Q 2 ) of 
the three terms in the integration X can approximately take their values at the same scale 
of Q 2 = k\. That is, the integration X is then simplified to 

1 1 3k 2 3 

1 = Jji-^M 1 + ak))fg(x] k 2 T ) = —-^.f g [x\ k 2 T ) = jpfgix; k 2 T ). (53) 



III. NUMERICAL RESULTS 

In this section we study the numerical behavior of the diffractive gluon jet production 
at the Fermilab Tevatron. We will study the pt distribution and x\ distribution of the 
cross section. We will also compare the gluon jet production with the quark jet production 
which has been calculated in [^,|T^]. A more thorough phenomenological study, including a 
comparison to currently available data at Tevatron on the diffractive dijet production rate, 
will be presented elsewhere. 
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Provided with the cross section formulas for the partonic processes qp — > qgp ( p23|) and 
ffP ~~ ► fl'S'] 5 (ID) we can calculate the cross section of the diffractive gluon jet production at 
hadron level. However, as mentioned above, there exists nonfactorization effect caused by 
the spectator interactions in the hard diffractive processes in hadron collisions. Here, we 
use a suppression factor to describe this nonfactorization effect in the hard diffractive 



processes at hadron colliders pjl0|. At the Tevatron, the value of J-'s may be as small 



as J-'s ~ 0.1 That is to say, the total cross section of the diffractive processes at the 

Tevatron may be reduced down by an order of magnitude due to this nonfactorization effect. 
In the following numerical calculations, we adopt this suppression factor value to evaluate 
the diffractive production rate. 

In our calculations, the scales for the parton distribution functions and the running 
coupling constant are both set to be Q 2 = k\. For the parton distribution functions, we 



choose the GRV NLO set 22 



In Fig. 4, we plot the differential cross section do j dt\ t= o as a function of the lower bound of 
the transverse momentum of the gluon jet, A/r min . This figure shows that the cross section is 
sensitive to the transverse momentum cut fcrmm- We plot separately the contributions from 
the two subprocesses, qp —>■ qgp and gp — ► ggp. By comparison, we also plot the cross section 
of the diffractive light quark jet production calculated in [TJJ]. The three curves in this figure 
show that the contribution from the subprocess gp — > ggp is two orders of magnitude larger 
than that from the subprocess qp — > qgp for the diffractive gluon jet production, and the 
light quark jet production rate is in the same order with that of the subprocess qp — > qgp. 
This indicates that the diffractive dijet production at hadron colliders dominantly comes 
from the subprocess gp — > ggp in the two-gluon exchange model. 

In Fig. 5, we plot the differential cross section do/dt\ t=0 as a function of the lower bound 
of the momentum fraction of the proton carried by the incident gluon Xi m i n , where we set 
&Tmin = 5 GeV . Fig. 5 (a) is for the contribution from the subprocess qp — > qgp, and Fig. 5(b) 
is from the subprocess gp —>■ ggp. These two figures show that the dominant contribution 
comes from the region of x\ ~ 10~ 2 — 10 _1 for the subprocess gp — > ggp, and x\ > 10 _1 for 
the subprocess qp — > qgp. These properties are similar to those of the diffractive charm jet 
and W boson productions calculated in fi~3|,[T5|j. 



IV. CONCLUSIONS 

In this paper, we have calculated the diffractive gluon jet production at hadron colliders 
in perturbative QCD by using the two-gluon exchange model. We find that the produc- 
tion cross section is related to the squared of the differential gluon distribution function 
dG(x; Q 2 ) I dlnQ 2 at the scale of Q 2 ~ k\, where kx is the transverse momentum of the final 
state gluon jet. We have also compared the production rate of the gluon jet in the diffractive 
processes with those of the light quark jet and heavy quark jet productions, and found that 
the production rates of these processes are in the same order of magnitude. 

As we know, the large transverse momentum dijet production in the diffractive processes 
at hadron colliders is important to study the diffractive mechanism and the nature of the 
Pomeron. The CDF collaboration at the Fermilab Tevatron have reported some results on 
this process 0. Up to now, we have calculated all of the dijet production subprocesses in 
the diffractive processes at hadron colliders, including gp — > qqp, qp —> qgp and gp — ► ggp 
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processes. In a forthcoming paper, we will compare the available data on the diffractive 
dijet production cross section at the Tevatron || to the predictions of our model to test the 
validity of perturbative QCD description of the diffractive processes at hadron colliders. 
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Figure Captions 



Fig.l. Sketch diagram for the diffractive dijet production at hadron colliders in perturbative 
QCD. The final state jet lines represent the outgoing quark or gluon lines, and the incident 
parton from the upper proton (labeled by Xip±) can be quark or gluon correspondingly. 
Fig. 2. The lowest order perturbative QCD diagrams for partonic process qp — ■> qgp. 
Fig. 3. The lowest order perturbative QCD diagrams for partonic process gp — > ggp. 
Fig. 4. The differential cross section da/dt\ t=0 for the gluon jet production in the diffractive 
processes as a function of £>rmin at the Fermilab Tevatron, where A>rmm is the lower bound 
of the transverse momentum of the out going gluon jet. 

Fig. 5. The differential cross section da/dt\t=o for the gluon jet production as a function of 
^imin, where Xi m i n is the lower bound of x\ in the integration of the cross section, (a) is for 
the contribution from the subprocess qp — > qgp, and (b) is from the subprocess gp — > ggp. 



17 




Fig.l 




Fig.2 



k + q 





Fig.3 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




I i I i I i I i I i I i I i I i I i I i I 

2 4 6 8 10 12 14 16 18 20 22 

Fig.4 




Fig.5(a) 




Fig.5(b) 



